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^ Abstract 

We investigate a method for determining the CP nature of a neutral Higgs boson or spin-zero 
resonance 4> at a future linear collider (ILC) in its $ — > t~t + decay channel. Our 
procedure is applicable if the production vertex of the Higgs boson can be measured. This 
will be the case, for example, for the Higgs-strahlung process e + e~ — > Z + 4>. We show 
that the method is feasible for both the leptonic and the hadronic 1 -prong tau decay modes, 
x ± ->■ Z ± + v T + V/, ->■ af, p ± , n ± -> n ± +X. 
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I. INTRODUCTION 



Recently, the ATLAS and CMS experiments reported the discovery of a neutral boson of 
mass ~ 126 GeV at the LHC 111 Q]. The experimental findings disfavor the option of a 
spin 7=1 resonance. The experimental results are compatible with the hypothesis of 
identifying this resonance with the Standard Model (SM) Higgs boson; however, much more 
detailed investigations will be necessary to establish this conjecture. The investigations of the 
properties of this resonance will probably be possible at the LHC to a large extent. 
A high-energy linear e + e~ collider would be an ideal machine to investigate the properties 
of this resonance, i.e., its couplings, decay modes, spin, and CP parity, in great detail (and, 
of course, also of other, not too heavy resonances of similar type if they exist). As it is likely 
that the ATLAS and CMS resonance is a spin-zero (Higgs) boson, one may revert, for as- 
sessing the prospects of exploring this particle at a future linear collider, to the many existing 
phenomenological investigations, within the SM and many of its extensions, of Higgs-boson 
production and decay in e + e~ collisions. As to the prospects of exploring the spin and CP 
properties of a Higgs boson, there have been a number of proposals and studies, including 



■ 2811 that are relevant for Higgs-boson production and decay at a linear collider. 
In this workshop contribution we apply a method [26, [27] for the determination of the CP 
properties of a neutral spin-zero (Higgs) boson $ in its T + T~ decays to the production of 
4> at a future e + e~ linear collider (ILC). For definiteness, we consider e + e~ — > Z4>, but the 
analysis outlined below is applicable to any other 4> production mode. In our analysis all 
major 1 -prong x decays are taken into account. We demonstrate that the CP properties of <I> 
can be determined with our method in an unambiguous way. 



II. CROSS SECTION AND OBSERVABLES 

Here we consider the production of a neutral Higgs boson 4» or, more general, of a spin- 
zero resonance of arbitrary CP nature by the Higgs-strahlung process in high energy e + e~ 
collisions: 

e + e~^Z + <t>. (1) 

For definiteness, we use = 126 GeV in the following. The following remark is in order 
here. As is well known, for a pure pseudoscalar boson $ = A, the AZZ vertex must be loop- 
induced 4 . We assume here, for the sake of choosing a definite production mode, that (OQ) 
applies also to the production of a pure pseudoscalar. 

For Z boson decays into an electron or a muon pair, the precise reconstruction of the produc- 
tion vertex and of the 4-momentum of the Z boson will be possible. As to <3>, we consider 



The strength of the loop-induced AZZ vertex was investigated for a number of SM extensions in 129 1. 
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here the decay mode into tau pairs, with subsequent 1 -prong x ± decays: 

t~t + ^a~a + +X, (2) 

where a = {e ,[l , K^} andX denotes neutrinos and, possibly, neutral pions. We assume 
that the tau-decay mode of the $ has a reasonably large branching fraction, which is the case 
in the Standard Model and in many of its extensions. The interaction of a Higgs boson <I> of 
arbitrary CP nature (J PC = ++ , J PC = h , or CP mixture) to x leptons is described by the 
general Yukawa Lagrangian 

& Y = ~(V2G F ) 1/2 m T (a T xx + b x xiy 5 x) <$> , (3) 

where Gf denotes the Fermi constant and a T , b T are the reduced x Yukawa coupling constants. 
In the following we take into account in © the main 1 -prong x decay channels 

x -> Z + V/ + v T , 

x -)> a\ +v T ->• n + 2K° + v T , 

x -> p + v T -> + tt + v T , 

x y k -\- v-c . (4) 

Our method that will be applied in the following does not require the knowledge of the x 
rest frame. Therefore we can include also the leptonic X decays in our analysis for which the 
presence of two, respectively four neutrinos would preclude the reconstruction of the x ± rest 
frames. We do not consider here x decays into 3 prongs, for instance x — > a\ — > 3 charged 
pions, because in this case the reconstruction of the x four-momentum should always be 
possible. This would considerably facilitate the measurement of the tau spin correlations that 
will be discussed below. (A corresponding analysis for production at the LHC was made in 
[Hi .) As an aside, we remark that it will be helpful, but not essential for future experimental 
analyses if the different hadronic T-decays can be experimentally distinguished. 
Our method to determine the CP properties of a spin-zero boson was first developed for the 



case of $ production in pp collisions at the Large Hadron Collider in 12611 and was then 
applied to an analysis that included the combination of all 1 -prong x decay channels in Ii27ll . 
The method is based on the fact that the CP quantum number of a neutral spin-zero resonance 
^> can be determined in a definite way through its x~x + mode by measuring the two x spin 
correlations S = s T - ■ s T + and Sep = k T • (s T - xs t +), where k T = k T /|k T | is the normalized 
T~ momentum vector in the zero-momentum frame (ZMF) of the T~T + -pair For a 

scalar the expectation value of S is (S) = 1/4, whereas for a pseudoscalar, (5) = —3/4. 
The CP-odd and T-odd spin correlation Sep probes whether or not 4> is a mixed CP state. If 
<J> is a CP mixture, i.e., if the neutral Higgs sector violates CP (that is, a T £<r ^ in ©), then 
a non-zero expectation value of Sep is generated already at tree level and can be as large 
as 0.5. 
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The spin of the x can not be measured directly; however it induces, in the spectrum of po- 
larized tau decay x ± — > a ± , a correlation with the direction of flight of the charged particle 



T(tT a T +x) dE aT dn a T/(47r) 

Here, s T denote the normalized spin vectors of the T T and q T the direction of flight of 
in the respective X rest frame. The function b(E a ) encodes the T-spin analyzing power of 
particle a. The correlation of the T-spins, s + ■ s~, leads to a nontrivial distribution of the 
opening angle Z(q + , q~), whereas the CP-odd observable k- (s + x s~) induces the triple 
correlation k • (q + xq"). The strength of these correlations depends, for a given strength of 
the reduced Yukawa couplings a Tl b T , on the product b(E a i-)b(E a +), while n(E a i-)n(E a +) is 
responsible for the rate of t + t~ decay into a + a'~ final states. 

In order to use these observables in an experimental analysis, one has to be able to reconstruct 
the x ± and a ± momenta in the x ± and 4> rest frames. This is, in general, not possible for the 
leptonic T-decay channel and very difficult in the case of hadronic x decays, because at a 
linear collider beamstrahlung effects can shift the initial center of mass energy by a large 
amount. In [26] it was shown that one can, nevertheless, construct experimentally accessible 
observables that have a high sensitivity to the CP quantum numbers of the The crucial 
point is to employ the zero-momentum frame of the a + a'~ pair. The distribution of the angle 

<p* = arccos(n^_ + • n"J_ - ) (6) 

discriminates between CP = ±1 states. Here are normalized impact parameter vectors 
defined in the zero-momentum frame of the a + a'~ pair. These vectors can be reconstructed 
Il26ll from the impact parameter vectors n T measured in the laboratory frame by boosting 



the 4-vectors n 1 ^ = (0,n=p) into the a'~a + ZMF and decomposing the spatial part of the 
resulting 4-vectors into their components parallel and perpendicular to the respective %^ 
or / T momentum. We emphasize that (p* defined in Eq. © is not the true angle between 
the x decay planes, but nevertheless, it carries enough information to discriminate between 
CP-even and CP-odd Higgs bosons. The role of the CP-odd and P-odd triple correlation 
introduced above is taken over by the triple correlation = p*L ■ (n^ + x n^~) between the 
impact parameter vectors just defined and the normalized a'~ momentum in the a'~a + ZMF, 



which is denoted by p!_. Equivalently, one can determine the distribution of the angle I126h 



V/^ = arccos(r-(Al + xnr)). (7) 

Before presenting results we would like to point out the difference of our method as compared 
to a previous analysis of how to determine the CP parity of a Higgs boson in its x + x~ decays 



at a linear collider. Refs. |20|-|22j] analyzed the hadronic 1 -prong decay x —¥ pv. The observ- 
able used by these authors, namely the acoplanarity angle of the p + and p decay planes, 
requires the reconstruction of the p + p _ ZMF, i.e., the measurement of the %^ and the n° 
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momenta, and the reconstruction of approximate x ± rest frames. As emphasized above, our 
method is applicable to all 1-prong x decays, in particular x — > I. 

III. RESULTS 

For predicting the distributions of the observables <p* and \jf^ P , for a specific 4>-decay mode 
©, in terms of the unknown reduced Yukawa couplings a t ,b t , one needs to know the spectral 
function n(E a ) and in particular b(E a ), which determines, as mentioned, the tau-spin analyz- 
ing power of particle a ± = / ± , 7i ± and therefore the shapes of the <p* and distributions. 
For the purpose of our analysis, the major 1-prong x decays © can be considered to be ex- 
perimentally well established Standard Model physics, and the respective spectral functions 



are known within the SM to sufficient accuracy, cf. |125|, 12711 and references therein. 
At a linear e + e~ collider, a Higgs boson produced by the Higgs-strahlung process (OQ) will 
have on average a much larger transverse momentum as compared to <J> production at the 
LHC by its maj or p roduction mode gg — > 4>. This calls for a study - independent of the 



LHC analyses II26L 12711 - of the question how the CP properties of a (pseudo)scalar boson 
are reflected in the distributions of (p* and i/z^p. In addition, for future experimental analyses, 
differences between LHC and ILC can be expected from the fact that the ILC detectors will 
be able to measure the x decay products at transverse momenta as small as about 10 GeV. 
As outlined above, our method is based on the reconstruction of the normalized spatial impact 
parameters of the T-decay products with respect to the production vertex of the Higgs boson. 
In the Higgs-strahlung process e + e~ — > 4> + Z at the ILC the normalized impact parameters 
can be reconstructed for events where the Z boson decays into electron or muon pairs, Z — > 
e + e~, and for <I> — > xx decays with sufficiently long T-decay lengths. Here, we use 

this process to study the (p* and \j/^ p distributions. We apply the following acceptance cuts 
appropriate for the ILC: 

pf > 10 GeV, 
15° < d l)7t < 165°. 

Let us first consider the decays x~x + — > l~l' + 4v, where = The functions n(Ej) and 
b(E[), where E\ is the energy of / in the x rest frame, are shown in Fig. The function 
b(Ej), which determines the correlation of the x spin with the momentum of /, changes sign 
at Ej = m T /4. For a Higgs boson with specified CP parity (and specified reduced Yukawa 
couplings), the sign of the product b(Ei)b(E t i) determines the functional form of the <p* 
distribution and in particular the sign of the associated asymmetry 

_ N(cp*>n/2)-N(cp*<n/2) 
r N((p* > n/2)+N((p* < k/2) ' K) 

As an illustration we apply the cut E[ > m T /4 in the x rest frames to allow only for contri- 
butions with b(Ej) < 0. The resulting normalized <p* distributions are shown as black dotted 
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Figure 1: (a) The spectral functions n(Ei) and b{E{), Eq. (f5]), for the leptonic T decays. The function 
n(Ei) is given in units of GeV -1 . (b) Normalized <p* distribution for a final state for a Higgs 
boson with mass of 126 GeV produced at yfs = 250 GeV. Scalar (j) = H, pseudoscalar (j) = A. The 
dashed lines show the distribution if no cuts are applied, the solid lines refer to the case where a cut 
p*-rest y 20 GeV was applied on both lepton energies in the Higgs rest frame. The dotted lines show 
the results for the ideal cut Ef~ iest > m T /4. 

lines in Fig.Q}). For a CP-odd boson it has its maximum at (p* = and its minimum at (p* = %, 
while for a CP-even boson (red dotted line) the distribution is flipped, (p* -h- % — <p*. If one 
applies instead, either for x + or T~ decay (but not for both), the cut E\ < m x /A which leads to 
b(E[) > 0, the behavior of the H and A distributions with increasing <p* will be interchanged. 
The magnitude of the resulting asymmetry © becomes smaller because the maximum of 
\b(Ei) \ is smaller for E\ < m T /4 than the maximum of \b(Ei) \ for E\ > m T /4. In addition, 
with the cut Ej < m T /4 the total decay rate is smaller than for E\ > m T /4. This would make 
a measurement more difficult. Without a cut on the lepton energy, the asymmetry of the 
normalized <p* -distributions is reduced, but remains non-zero because the averaging over the 
different signs of b(Ei) is weighted by the spectral function n{Ej) displayed in Fig. [T^. The 
result is shown in Fig. [TJ? (dashed lines, black for a pseudoscalar, red for a scalar boson). 
Obviously, it would be an advantage if one could apply a cut on Ej to separate the contribu- 
tions that involve different signs of b{E{). However, this would require to reconstruct the full 
T 4-momentum in order to perform the necessary boost into the T rest frame. On the other 
hand, the energy of the lepton / in the Higgs rest frame, Ef, is correlated with the energy E\ 
in the T rest frame and a cut on the former can enrich the event sample with events in the de- 
sired range of the latter. The Higgs rest frame can, in fact, be reconstructed for the production 
process e + e~ — > <t> + Z, because the Z-boson 4-momentum is known for Z — > e, /i decays, 
provided that initial state radiation is negligible or can be corrected for. This should be the 
case at least for the TESLA design [30]. In Fig. [Th we show the resulting <p* distributions 
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with a cut Ef > 20 GeV applied to both leptons from T decay. The solid black line (solid 
red line) shows the distribution for a CP-odd (CP-even) Higgs boson. The sensitivity of the 
distributions to the CP parity of <E> is clearly enhanced compared to the case where no cut is 
applied (dashed lines). These distributions are only slightly less sensitive than those with the 
ideal cut Ef~ KSt > m T /4 (dotted lines). 

In the 2-body decay x — y Tt -\- Vf, the % is monochromatic in the x rest frame. (Its T-spin 
analyzing power is maximal.) A cut on the energy of the charged prong, i.e., of the charged 
pion, is very important for the hadronic 1 -prong decays x — > p ± v — > iz^ + 7T°v and x ± — > 
afv — > Ti^ + 2k°v. For example, for the decay p^ — > 7T =t + 7t°, the function b(E K ) changes 
sign within the range of E n (see, e.g., Fig. 4a in fl270) and the <p* distributions for both a scalar 
and a pseudoscalar boson turn out to be flat if no cut was applied. The same is true for the 
x — > a\ decay mode. As in the case of x — > I, a cut on the energy of the charged pion in the 
Higgs-boson rest frame such that b(E n ) is either positive or negative for the selected events 
significantly enhances the discriminating power of the (p* distribution. Provided that the event 



rate is large enough, a value for E^ cut may be chosen such as to optimize the separation of 
events with positive and negative b{E n ), and both sets of events could be used to determine 
the CP nature of <J>. 
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Figure 2: Influence of a cut on the energy of the charged pion in / + p final states for a scalar boson 
(j> = H with a mass of 126 GeV produced at y/s = 250 GeV: (a) on the <p* asymmetry, (b) on the cross 
section. The black full lines correspond to applying the cut Ef < Ef cut , while the red dashed lines are 
for E° > 

We illustrate this in Fig. [2]for the decay 4> — > I p + + 3v of a scalar boson. The effect can be 
quantified by calculating the associated asymmetry ©. The dashed red line in Fig. [2^ shows 
the asymmetry for events with > £* cut . The value Ef cut = 0, not shown in the figure, 
corresponds to the case without cut. In this case the asymmetry is rather small; applying 
a cut, the asymmetry increases to almost —12% for Ef cut = 35 GeV. However, increasing 
£*cut w iU decrease the cross section as shown in Fig. [2}). (The cross section was computed 
for the Standard Model Higgs boson.) Without any cut the cross section is about 1.5 fb; it 
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Figure 3: Normalized \j/^ p distribution for Ip final states for different types of Higgs bosons with 
mass of 126 GeV produced at y/s = 250 GeV. A minimum cut on the 7T + energy of 25 GeV in the 
Higgs rest frame was applied. 



decreases to 0.2 fb for the cut Ef cut = 45 GeV. 



On the other hand, the complementary region < E^ cut leads to a positive asymmetry Ay* 
for cut values < 35 GeV due to the fact that b(E n ) has changed sign. For small values of E n , 
one finds A<p* ~ +9%, but the cross section is tiny, about 0.05 fb. The asymmetry decreases 
to almost —2% for Ef cut = 45 GeV. It is clear that a judicious choice of this cut, taking into 
account experimental conditions and the available luminosity, is required to reach an optimal 
discrimination between a CP-even and CP-odd boson. 

The (p* distribution is well suited to distinguish between CP-even and CP-odd states. How- 
ever, if the Higgs boson is a CP-mixture, or if there would exist two (almost) mass-degenerate 
bosons that escape experimental resolution, one of which has CP parity +1 and the other one 
— 1, the (p* distribution would be flat, assuming the cross sections are of comparable magni- 
tude. The distribution with respect to the angle defined in Eq. © would be appropriate 
to resolve these scenarios. A typical result is shown in Fig.[3]which applies to the decay chain 
4> — > T + T~ — > l~K + + 3v via hadronic T — > p decay. The solid blue line shows the normal- 
ized \y* distribution of a maximally CP-mixed boson (\a T \ = \b T \ > 0) and we have chosen 
a r = —b T . For the scenario of mass degenerate bosons with opposite CP parities the 
distribution is shown by the horizontal dashed black line. The dotted blue line corresponds to 
the case of a non-maximal mixture with reduced Yukawa couplings a x = —2b T . If a T and b T 
have the same sign, the distribution is flipped, \j/^ p -h- % — ^q P - The resulting asymmetry of 
the \\Tqp distribution will clearly be observable, provided the event rates are large enough. 
We have also performed a Monte Carlo study to estimate the uncertainty of the measurement 
of the impact parameters. We applied a simple Gaussian smearing with compact = 

25° \m 

on the direction of the normalized impact parameter vectors. We found that the <p* and 
distributions are only mildly affected by such an uncertainty. 
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IV. CONCLUSIONS 



Using the method of fl26L 12711 we have shown that the CP nature of a neutral Higgs boson $ 
produced at a future linear collider can be determined in a definite way in the <I> — > T + T~ 
decay channel with subsequent 1 -prong x decays. We have considered the production of <E> 
with mass m<j> = 126 GeV by the Higgs-strahlung process, but our method can also be applied 
to Higgs bosons with other masses and to any other $ production mode. Our approach does 
not require the knowledge of the x rest frames; therefore, all 1 -prong x decays, including 
x — > I can be used. The joint measurement of the distributions and associated asymmetries of 
the angles (p* and y/£p allows to discriminate between a number of scenarios, some of which 
were discussed above. For a statistically significant determination of the CP parity of 4>, only 
very few 1 -prong events are required. We will elaborate on this and on other related issues in 



future work [311. 
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